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Abstract A knot-theoretic explanation is given for the rationality of the quenched QED 
beta function. At the link level, the Ward identity entails cancellation of sub divergences 
generated by one term of the skein relation, which in turn implies cancellation of knots 
generated by the other term. In consequence, each bare three-loop diagram has a rational 
Laurent expansion in the Landau gauge, as is verified by explicit computation. Compa- 
rable simplification is found to occur in scalar electrodynamics, when computed in the 
Duffin-Kemmer-Petiau formalism. 
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1. Introduction The surprising rationality of the three- |rj and four- loop 
quenched (i.e. single-electron-loop) terms in the QED beta function is an outstanding 
puzzle, which we here elucidate by giving the Ward identity Z\ = Z2 an interpretation in 
terms of the skeining relation that is the basis of the recent association of knots |3], |], |5], |J 
with transcendental counterterms. 

In Section 2, we study the intricate cancellations between transcendentals in all 6 
of the methods (known to us) for calculating /3(a) = da/dln/1 2 = J2n Pnd n+1 , with a 
QED coupling a = a /Arc. Our focus is the rationality of = — 2 || and f3± = 
—46 0. (Subscripts denote the number of loops; where necessary, superscripts denote 
the number of electron loops.) In Section 3, we compound the puzzle by exposing even 
more intricate cancellations of £3 in three-loop scalar electrodynamics 0]. The argument 
from knot theory is given in Section 4, leading to specific predictions, confirmed by detailed 
calculations, performed using the techniques of || ||, [UJ. Conclusions, regarding higher 
orders 0, unquenched (i.e. multi-electron- loop) contributions [|TT], |TJ, [T3|] , and non-abelian 
gauge theories O, |j~5|, are presented in Section 5. 



2. Six calculations in search of an argument There are (at least) 6 
ways of obtaining the beta function of quenched QED. For each, we expose the delicate 
cancellation of transcendentals between diagrams, which cries out for explanation. 

Method 1: Dyson-Schwinger skeleton expansion |l|, 16| . The Dyson-Schwinger equa- 
tions give the photon self-energy, schematically, as |T7| 

n> = ZjT^ = ry?(i - kg)t u , (1) 

where T M is the dressed vertex and G stands for the pair of dressed propagators in Fig. la. 
To obtain the second form, illustrated in Figs. lb,c, one uses T u = Zi^ v + KGT U , where K 
is the kernel for e + e~ scattering and loop integrations and spin sums are to be understood 
in the products. Now we expand each vertex to first order in the external momentum q: 
r /t = T° + q ■ dr^ + 0(q 2 ) where d a = d/dq a and the Ward identity gives T° = T^(p,p) = 
{d/dp fJ- )S^ 1 (p) in term of the inverse propagator. We use the low-momentum expansion 

(r„ - rJ)G(i - kg)(t u - r°) = ( q ■ dr„)G(i - KCf)( q ■ dv w ) + o( q 3 ) , (2) 

differentiate twice, and make liberal use of (1 — KG)d a Y v = d a (KG)T® + 0(q), to obtain 

'ld a dfsG + (d a G)K(d^G) + 2{d a G){df]K)G + \G{d a d p K)G 
+ (d a (GK))G(l - KG)-\dp(KG)) + (a ^ /3)1 T° u + 0(q) , (3) 



d^IV = 



which entails only K and S. There is a very simple statement of this result: between 
the dressed zero-momentum vertices occur all and only the terms in d a d{3(G(l — KG)^ 1 ) 
that give no subdivergences. In other words, the Dyson-Schwinger method kills maximal 
forests of subdivergences both on the left and on the right; only the overall divergence 
survives. Moreover, after cutting at a line on the left (or right), we may set both the 
mass and the external momentum to zero, with no danger of infrared divergence, and 
hence obtain the L-loop quenched beta function from finite massless two-point skeleton 
diagrams, with up to (L — l)-loops. 

We route half of the external momentum through the electron line and half through 
the positron line and find that each of the 5 terms in (|3]) yields a contraction-independent 



1 



contribution to /3 3 
where 7 = 7 2 °' = 



[i] 



To all orders, the first contribution to /3^(a) is |a 2 (l 



|7 2 + 7 3 ) 



£a 



3^2 



+ |a 3 + 0(a 4 ), is the quenched electron-field anomalous 



dimension |18j with a photon propagator g^ u /k 2 + (£ — ^k^k^/k 4 . The second term gives 
/3 2 = 4. At 3 loops, all 5 terms contribute, giving a remarkable cancellation of £3: 



4 1] 



-2£ 2 



+ 
+ 



(7£-3)(£-3) + 2(£ 2 -6£-3)C 2 



j(£ 2 + 12£-23) + 2(£ 2 -6£ + 13)C2 



+ 8 = -2. 



5) f-4( 3 



(4) 



In the Landau gauge, £ = 0, we reproduce Rosner's (^-cancellation: = — 6 — 20 + 26 |T[. 
In the Feynman gauge, £ = 1, the third term vanishes; the second and fourth are still 
transcendental and exactly cancel the fifth. 

Method 2: Integration by parts of massive bubble diagrams ||. We may study sepa- 
rately each of the 3-loop quenched diagrams of Fig. 2, using the techniques of ||, which 
reduce to algebra the calculation of 3-loop massive bubble diagrams, in d = 4 — 2e 
dimensions, thereby also yielding the finite parts of on-shell charge renormalization, 
used in ]I0| to establish the connection between 4-loop on-shell and minimally sub- 
tracted (unquenched) beta functions. The coefficient, d n , of a n /e, in the sum of all 
n-loop, single-electron-loop, bare diagrams contributing to 1/Z 3 , differs from (3$/n; one 
must also take account of the quenched anomalous mass dimension, 7^ = I] n 7 n a n = 
3a + fa 2 + ^ a 3 + 0(a 4 ) @. To 4 loops, the effects of mass renormalization are also 
rational, giving agreement with 3-loop results || [TQfl : 



di = Pi 



|/3 2 -2/3 l7l 



-6; 4 = I/?W-2/3 2 7i-/5 1 72 + 6/? 1 7 1 2 = ^; (5) 



and generating a 4-digit prime in the quenched 4-loop coefficient: 



d A 



&72 + mil + 8/3i727i - IA73 - f A7i 3 



2969 
~6~ 



(6) 



which string-inspired techniques [19] may eventually reproduce. In the simplest case || of 
contracting dadpU^lq^ with g^ga/3 in the Feynman gauge, we find cancelling coefficients 
of Cs/e in the contributions of Figs. 2f and 2g, corresponding to (^-cancellation between 
the fourth and second terms, respectively, of (f|) at £ = 1. 

Method 3: Integration by parts of massless two-point diagrams p0| , 
study the behaviour of the bare diagrams of Fig. 2 at large q 2 
m = and have no need of mass renormalization, or differentiation. The integration-by- 



Next we 



21], where we may set 



parts method of Chetyrkin and Tkachov []20| now suffices. It has been implemented in the 
program mincer, whose test suite [32"| evaluates all the diagrams of Fig. 2 in the Feynman 
gauge. Table 1 of |22| reveals that, in addition to the crossed-photon diagrams of Fig. 2f 
and 2g, the uncrossed-photon diagrams of Figs. 2a and 2e also entail (s/e. Moreover, the 
cancellation does not occur in the crossed and uncrossed sectors separately; the relative 
weights of (3/e in the contributions of the bare diagrams 2a,e,f,g are 3 : —6 : —1:4. 

Method 4 : Infrared rearrangement of massless diagrams |l4l |23], 24 1. Infrared rear- 



rangement is a technique for reducing the calculation of L-loop counterterms to that of 
(L — l)-loop massless two-point diagrams, by the subtraction of sub divergences in the MS 
scheme, followed by nullification of external momenta and appropriate cutting of massless 



2 



bubble diagrams. As shown by Method 1, the Dyson-Schwinger equations make this un- 
necessary in QED. However, the technique prospered at the 3-loop level in QCD [|TJ], ^3 
to such an extent as to encourage the calculation of the 4-loop beta function of QED, with 
the rational quenched result /3| = —46 . A measure of the seemingly miraculous cancel- 
lation of transcendentals in this method is afforded by the complicated combination of £3 
and (5 that occurs in the non-abelian 4-loop QCD corrections to R(e + e~ — *> hadrons) 



Method 5: Propagation in a background field p6|, p7|. As observed in |2£j, the derivative 
(d/ da) (a) / a 2 ) = Y, n >i( n ~ l)/^ ™ -2 ma Y be obtained from the coefficient of aF^ u in 
the single-electron-loop contributions to the large-g 2 photon propagator in a background 
field Fnv, thereby providing a further alternative to infrared rearrangement. Two-loop 



massless background-field calculations have been performed in QCD [E7|, with results 
recently confirmed by a full analysis of the massive case p9[. From the two-loop cor- 



rections, 1 + (2Ca — Cp)a s /47r, to the coefficient of <a s G > in the correlator of the 



light-quark vector current of QCD, one immediately obtains 2(3^' /(3 2 = —1, by setting 
Ca — 0, Cp — 1, which confirms the correctness of this method at 3 loops. A measure 
of the complexity of the (3 cancellations is afforded by studying the diagram-by-diagram 
analysis of the appendix of p8 . 



Method 6: Crewther connection to deep-inelastic processes |30, jflfl . Finally, a useful 



checkoff 11 = -46 was obtained in |5D| by taking the reciprocal of the 3-loop |JT| radiative 
corrections to the Gross-Llewellyn-Smith sum rule, in the quenched abelian case: 







Mi 



4 2 



1 



3a + fa 2 -fa 3 + 0(a 4 ) 



fa 2 + 4a 3 - 2a 4 - 46a 5 



0(a b 



(7) 



in precise agreement with 0. This is an example of a Crewther connection |32fl which 
is unmodified by renormalization in quenched QED. A measure of the complexity of the 
cancellations of ( 3 and ( 5 is afforded by studying the non-abelian terms in 3-loop deep- 
inelastic radiative corrections which are replete with both transcendentals. 

3. Cancellation of (3 in scalar electrodynamics Corresponding cal- 
culations in scalar electrodynamics (SED) are most conveniently performed using the 
Duffin-Kemmer-Petiau [|33[ (DKP) spinorial formalism for the charged scalar field, in 
which the Feynman rules are identical to those of fermionic QED, with a bare vertex eo7 M 
and a bare propagator S (p) = \j{j> — m ). The only difference resides in the 7-matrices: 
in the DKP formalism #i is not invertible; instead one uses the fact that = p 2 p* to obtain 
Sq(p) = {^{i> + m o)/(p 2 ~ m o) ~ l)/ m o- The trace of the unit matrix is d+ 1, the trace 
of an odd number of 7-matrices vanishes, and the trace of an even number, 7 /Ltl to 7 M2n , 
is the sum of two terms [j3"4| : g MlAt2 • • • Q^-i^n an d the cyclic permutation a M2/tl3 . . . g^^- 
The effect is to make all traces regular as mo — > 0, while automatically generating the 
many seagull terms of conventional scalar methods [17|. This is a great simplification, 



eliminating the need to include 21 seagull diagrams at 3 loops, which would lead to many 
more terms in d a dpYi^ u \ q= Q. 

We compute the diagrams of Fig. 2 using the DKP formalism. At q = (Method 2) 

program RECURSOR ||, for 3-loop massive bubble diagrams; at 



we use the REDUCE 35 



m = 0, we use the reduce program slicer [10], devised to check the results of [21] for the 
large-g 2 photon propagator in the MS scheme. As in QED, we perform both calculations 
in an arbitrary gauge and contract d a d^II Ml/ with g^ u g a p + X(g^9a0 + 9^9vp + 9va9np), 
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where A is an arbitrary parameter, which affects the contributions of individual diagrams, 
but not the total result for the transverse self energy li^q) = (q^q u — q 2 g liV )Tl{q 2 ). At 
m = 0, the corresponding freedom is to contract with the tensor g^ u q 2 + X(d + 2)q^q u . 
On-shell SED mass renormalization is performed as in H [| : there are only two quenched 
self-energy diagrams at two loops; each is projected on-shell by taking the trace with 



+ m) at p 



m 



with a pole mass m. Unlike the QED case, the relation between 



bare and pole masses is infrared-singular in SED, though that causes no problem for the 
dimensionally regularized calculation of Z 3 = l/(l + n (0)), where infrared singularities in 
the bare diagrams for IIo(O) are cancelled by those in Z m = mo/m. The on-shell methods 
of !| then yield the quenched contributions 




e(2-e)(l-4e 2 ) 



-| + 16C 2 (5-81n2) + iC 3 + 



3c 
29 



+ f + 8C 2 (3-41n2) + f|C3 + 



77 
9 

136 
9 



for QED [^, [Tt| and SED, respectively, where a m = r(l-|-e)eg/ (47r) d//2 m 2<: is a dimensionless 
coupling, with a pole mass m, and terms of order af^e and are neglected. (There is no 
need to renormalize the bare charge, eo, when dealing with the quenched contributions.) 
From the singular terms in ([8]) we read off the 3- loop beta function of quenched SED: 
fiW (a) = \a 2 + 4a 3 + ^a 4 + 0(a 5 ). We have calculated the double-bubble term /3 3 2 ' in the 



on-shell and MS schemes, obtaining agreement with J7| 
For the 3-loop quenched MS contributions to TI(g 2 ) 



in the latter, as shown in Table 1. 
, at large q 2 , we obtain 



n 3 V 



-21n(/. 2 /g 2 )-^-^C3 + 160C 5 
^ - ^ C 3 + 40C 5 



29 
2 



ln(fi 2 /q 2 ) + 



(9) 



with a = r(l+e)r 2 (l-e)eg/r(l-2e)(47r) d / 2 /i 2e , which suppresses (2 in Laurent expansions. 



The QED result confirms [^TJ. The SED cancellations are even subtler: with £ = A = 
the relative weights of from bare diagrams 2a,e,f,g,h are 18 : —36 : —5 : 22 : 1 . 

4. The argument from knot theory From the point of view of knot theory, 
as proposed in 0, |], ]5], 0, the presence of transcendentals in counterterms can be traced 
to the knots that are obtained by skeining the link diagrams that encode the intertwining 
of loop momenta in Feynman diagrams. The absence of transcendentals in the quenched 
beta function does not, therefore, correspond to the absence of knots in the Feynman 
graphs, since the crossed-photon graphs of Figs. 2f,g,h all realize the link diagram whose 
skeining contains the trefoil knot ||. Accordingly we expect to find (3/e in their divergent 
parts. To explain the cancellation of transcendentals, we must study the interplay between 
knot-theoretic arguments and the gauge structure of QED. 

It was found in 0, |J that ladder topologies are free of transcendentals when the 
appropriate counterterms are added: after minimal subtraction of subdivergences, ladder 
graphs, such as in Figs. 2a,e, give rational terms in the Laurent expansion in powers of 
1/e. In [EL H], on the other hand, transcendentals corresponding to positive knots, with 
up to 11 crossings, were successfully associated with sub divergence-free graphs, up to 7 
loops. Thus the skein relation played two distinct roles in previous applications: in || 
the so-called A part of the skein operation determined the subdivergences, while the B 
part gave no non-trivial knots in ladder topologies; in || there were no subdivergences 
associated with the A operation, while the knots from the B operation faithfully revealed 
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the transcendentals resulting from nested subintegrations. In Figs. 2f,g,h we are now 
confronted with Feynman diagrams whose link diagrams generate the trefoil knot (via B) 
and also have sub divergences (corresponding to A). 

We thus propose to associate the cancellation of transcendentals with the cancella- 
tion of sub divergences in the quenched beta function of QED, which is an immediate 
consequence of the Ward identity, Z\ — Zi. 

To see the key role of the Ward identity, consider Fig. 2g. There is an internal vertex 
correction, which is rendered local by adding the appropriate counterterm graph. Due 
to the Ward identity, this counterterm graph is the same as that which compensates for 
the self-energy correction in Fig. 2e. In |3|, |J it was shown that the latter counterterm 
could be interpreted as the A part of the skein operation on the link diagram L(2e) of 
Fig. 3, associated with the Feynman graph 2e. We assume that this is a generic feature 
of the relationship between skeining and renormalization and associate the corresponding 
counterterm for Feynman graph 2g with the term obtained from applying A twice to 
the link diagram L{2g\ which requires two skeinings to generate the same counterterm, 
along with the trefoil knot from the B term. The Ward identity thus becomes a relation 
between crossed and uncrossed diagrams, after skeining: 



The trefoil knot results, in this language, from B(B(L(2g))), which generates, in general, 
a (3/e term from Fig. 2g, even after the subtraction of subdivergences. We now use the 
Ward identity ([H]) at the link level to obtain 



from which we see that it relates the transcendental counterterm from a torus-knot |4[] 
topology, in Fig. 2g, to a knot-free [[| ladder topology, in Fig. 2e. 

We conclude that should be absent from the bare diagram of Fig. 2g when it is 
calculated in the gauge where the bare diagram of Fig. 2e is free of C3A, i-e. in the Landau 
gauge, where the latter is free of subdivergences. In Fig. 3 we summarize the argument. 
By a similar argument, we also conclude that the other graphs with the trefoil topology, 
namely Figs. 2f,h, should be free of in the Landau gauge. 

Note that vanishing of the bare ladder diagram of Fig. 2e, in the Landau gauge, does 
not imply the vanishing of B (B (A^ 1 (L(2e)))); it merely implies its rationality. The action 
of the A and B operators is non-trivial on the diagram of Fig. 2e; in the language of 0, |j, 
skeining involves a change of writhe number in the sub divergence, and only terms with 
writhe number zero vanish in the Landau gauge, while the double application of the B 
operator in B(B(A~ l (L(2e)))) generates a non-vanishing writhe for the sub divergence. 

To test these ideas, we evaluate the separate contributions of the QED diagrams of 
Fig. 2 to (||,^D in the Landau gauge and find, indeed, that each bare diagram has a rational 
1/e term in its Laurent expansion, at £ = 0, for any value of the contraction parameter, 
A. In any other gauge, the bare massless diagrams of Figs 2a,e,f,g give non-zero, mutually 
cancelling, coefficients of Cs/ e - The m = contributions to the Laurent expansion of LT^ 
are given in Table 2, in the case £ = A = 0. The same rational behaviour is observed in 
the Landau gauge at q = (again for any value of A) with contributions also recorded in 
Table 2 (with A = 0), along with that from mass renormalization of lower- loop diagrams. 



A(A(L(2g))) = A(L(2e)) A(L(2g)) = L(2e) . 



(10) 



B(B(L(2g))) = B(B(A- 1 (L(2e)))), 



(11) 
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To determine the corresponding gauge for SED, we study both of the two-loop DKP 
diagrams and find that each gives 1/e 2 terms, proportional to 2£ — 3. Thus the gauge 
£ = 3/2 is the closest DKP analogy to the Landau gauge of QED. In any gauge, however, 
diagram 2h gives a (3/e term, proportional to A — 2. (The absence of such a term in QED 
is attributable to the fact that the finite part of the two-loop fermion propagator does 
not involve (3 at large momentum.) Accordingly we evaluate all m = and q = SED 
diagrams with £ = 3/2 and A = 2 and find that every one is indeed rational, with a Laurent 
expansion recorded in Table 3. Again, we regard this as a significant success of skeining 
arguments, since fixing 2 parameters removes (3/e from each of the 5 diagrams 2a,e,f,g,h. 

5. Conclusions Complete cancellation of transcendentals from the beta function, 
at every order, is to be expected only in quenched QED and quenched SED, where sub- 
divergences cancel between bare diagrams. 

Contributions with more than one charged loop require scheme- dependent coupling- 
constant renormalization. Knot theory guarantees that double-bubble contributions are 
rational at the three- loop level, as found in |7], 1 1]| and shown in Table 1. But at four 



loops 0, [n], [T2[], and beyond |T3|] , such multi-electron- loop contributions entail non-trivial 



knots, associated with zeta functions and more exotic |5|, |(J transcendentals, whose 
cancellation is not underwritten by the Ward identity. 

Likewise, we do not expect the beta functions of non-abelian gauge theories to be 



rational beyond the presently computed [13, [TH] three-loop level, since coupling-constant 



renormalization is required to remove sub divergences generated by gluon- and ghost-loops 
in the gluon self-energy, which cannot be quenched without violating the Slavnov- Taylor 
identities. Indeed, it is still somewhat obscure how the rationality of /?3 comes about in 
QCD: as shown in ||15| , the renormalization constants Zi.2,3, for the quark-gluon vertex, 
quark field, and gluon field, all involve (3/e at three-loops; only in Z a = ZfZ^Z^ 1 does 
the cancellation occur. It may be possible to understand this in a three-loop background- 
field calculation, but we do not expect the intrinsically scheme-dependent non-abelian 
beta function to remain rational at higher orders in the MS scheme. 

In fact, there are only two further quantities that we expect to be rational beyond three 
loops: the quenched anomalous field and mass dimensions, 72°' and 7^, of QED, whose 
rationality at the three-loop level is established in ITR 113. These are scheme-independent, 



since sub divergences cancel in the electron propagator, by virtue of the Ward identity, 
which here relates proper self-energy diagrams to one-particle-reducible diagrams. 

One is tempted to seek a calculational method, based on the skein relation, to compute 
the rational anomalous dimensions of quenched QED, in terms of diagrams, with non- 
zero writhe numbers, that have been 'ladderized' by Ward identities, such as (|ll|), which 
may eventually yield a rational calculus like that in [|J. The complexity of the rational 
contributions of Tables 2 and 3 indicates that such a calculus would be rather non-trivial. 

In conclusion: all-order rationality of counterterms is specific to quenched abelian 
theories, where the cancellation of knots, from the B term H |6| in the skein relation, 
matches the cancellation of sub divergences, from the A term || ||. This is exemplified 
by the rationality of all quenched three-loop bare QED diagrams, in the Landau gauge. 
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Table 1: Three- loop on-shell and MS coefficients of QED and SED beta functions. 





QED 


SED 


on-shell 


-2-S* = -3H [11] 


29 55 _ 103 

2 18 9 


MS 


-2-f = -f m 


29 49 _ 106 m\ 
2 18 9 HO 



Table 2: QED contributions to n 3 , at m = and g = 0, with £ = A = 0. 



method 


m = 


q = 


Laurent term 


a 3 /e 2 


a 3 /e 


</e 2 


a 3 Je 


diagram 2a 





8/3 





8/3 


diagram 2b 








-12/5 


152/5 


diagram 2c 








72/5 


8/5 


diagram 2d 








2/3 


31/3 


diagram 2e 











-24 


diagram 2f 


-2/3 


7/3 


-2/3 


17/3 


diagram 2g 





-28/3 





-28/3 


diagram 2h 


2/3 


11/3 


-12 


28 


uncrossed 





8/3 


38/3 


21 


crossed 





-10/3 


-38/3 


73/3 


mass renorm 











-46 


total 





-2/3 





-2/3 



Table 3: SED contributions to II3 , at m = and q = 0, with £ = 3/2 and A = 2. 



method 


m = 


q = 


Laurent term 


a 3 /e 3 


a 3 /e 2 


a 3 /e 


</e 3 


a 3 Je 2 




diagram 2a 


1/16 


55/96 


4151/576 


1/16 


7/96 


2501/576 


diagram 2b 


-13/8 


-409/24 


-16243/144 


-13/8 


-233/120 


-8/45 


diagram 2 c 


3/2 


31/2 


3703/36 


3/2 


43/20 


203/45 


diagram 2d 


1/16 


103/96 


4387/576 


1/16 


-65/96 


-131/576 


diagram 2e 





-5/48 


41/32 





-5/48 


-1141/96 


diagram 2f 


-1/16 


-9/32 


911/192 


-1/16 


7/32 


1013/192 


diagram 2g 





3/16 


-917/288 





3/16 


-1349/288 


diagram 2h 


1/16 


3/32 


-1667/576 


1/16 


-29/32 


-809/576 


uncrossed 








37/6 





-1/2 


-55/16 


crossed 








-4/3 





-1/2 


-13/16 


mass renorm 














1 


109/12 


total 








29/6 








29/6 
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Fig. 1: Illustration of terms in the Dyson-Schwinger equation (1). 

O O <H> 

a b c 



Fig 2: Diagrams contributing to /3 3 , with photon lines drawn inside the electron loop. 




e f g h 

Fig. 3: The A part of the skein operation on the link diagrams for Figs. 2e,g. 




2e L(2e) A(L(2e)) 




2g L(2g) A(L(2g)) = L(2e) 



